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Abstract 

We study an inflationary model where two decoupled string axions drive inflation. The 
number of e-folds is dependent on the energy scale and the decay constant, but is almost 
independent of the angular component in spite of the rich geometry of the moduli space. 
This suppresses the nearly scale-invariant spectrum of the isocurvature perturbations, 
making the power spectrum of the primordial density perturbations dominated by the 
adiabatic component. We also briefly discuss some amendments to improve the situation. 
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1 Introduction 



Inflation [lj is currently the most promising candidate to solve many cosmological problems. 
It is, however, still unclear yet how to implement a consistent and successful scenario in the 
context of string theory. Fortunately, after the de Sitter vacua construction in type IIB theory 
[2], there have been considerable advances to this end for recent years [3|. An inflaton candidate 
of particular interest is the string axior0, which is known to exist abundantly in string theory, 
and they have flat potentials even after all the moduli are stabilised. 

In this paper, we study a simple inflation model which consists of two axion fields. Since 
there are plenty of axions, it is not difficult to construct an inflation model with the number 
of e-folds far larger than 60 using many fields. As inflation proceeds, however, fewer and fewer 
fields are responsible for inflation because heavy fields drop out of the inflationary regime and 
contribute no more. Thus, it is very conceivable that in the final stage of inflation, i.e. the 
last 60 e-folds relevant for the observable universe, significantly smaller number of fields drive 
inflation. In this regard, it is natural to consider the simplest two-field model first which is 
considered important in building inflation models with multiple fields [5]. 

This paper is outlined as follows. In Section [2] we briefly review the potential of two axion 
fields. In Section [3] we study the dependence of the number of e-folds on the parameters of the 
potential and the density perturbations produced during inflation. In Section @] we discuss how 
the problems associated with the simplest model could be alleviated. In Section \5\ we give a 
brief summary. 



2 Potential 

As mentioned in the previous section, we will be interested in the two-field case, so the potential 
is written as 
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where and are the dynamically generated axion potential scales, and and are the 
axion decay constants. We expect both and f^) are smaller than mpi, especially the 

potential scale could be significantly so. There are many maxima, minima and saddle points, 
and they lead to a very rich structure of the potential depending on the values of the parameters. 
We illustrate a region of the potential in Fig. [T] 

For simplicity, we concentrate on = and A^, = A^, and occasionally we will omit the 
subscripts. Then from above we can write the Hubble parameter as 
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Near an extremum point, the effective mass is given by \m\ 

mpi 



2nA 2 /f. Hence 
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i.e. the mass is very heavy and therefore the slow-roll condition |m P1 V'/V| <C 1 is badly broken. 



1 See, e.g. Refs. [4] for detailed accounts on the cosmological aspects of axion. 
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Figure 1: (Left) plot of the potential, rescaled to V max = 1, and (right) contour plot of the same 
region. 

3 Inflation 

To compute the evolution of the scalar field on the potential given by Eq. ((TJ), it is convenient 
to modify the equations of motion for and if) and to eliminate the scale factor a(t). Using 
d/dt = Hd/dN, we can derive the equations in terms of the number of e-folds N as 

where a prime denotes a derivative with respect to N, Vj = dV/difii and the indices % and j 
denote the canonical fieldal and ip. 

Now we study the trajectory of the field starting from a maximum of the potential When 
the field is about to roll down the potential from a maximum, we need to specify two quantities, 
namely, the displacement (how far is the field off from the maximum?) and the misalignment 
(where is the field heading for?). It is plausible to expect that due to quantum fluctuations 
the field is off the maximum by an amplitude of H ~ A 2 /mp\. For the misalignment 9, which 

2 For non-canonical fields, i.e. 

where G^ is the metric of the field space, the equations are given by 

/' + 3/ - + ll^V' + ^ (6*& - G^ V) , 

where Tj k is the Christoffel symbol constructed by G^ . 

3 We don't consider the inflationary phase which occurs as the field rolls down from a saddle point. This 
saddle inflation is extensively studied in [B]. 
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should be written as 



t&nO = — 



(5) 



first we note that near a maximum we have an approximate U(l) symmetry. Hence we can infer 
that as long as there is enough time for the field to get arbitrary orientation, every direction is 
equivalent near a maximum. To estimate this time, we note that the field is placed away from 
the maximum by H. Therefore, we can change the question to ask how long it takes for the 
field to get randomly placed on a circle of radius H. To have inflation, we require that 



and hence 



Thus, we can estimate 
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i.e. the field is randomly placed almost instantaneously. 

Some results of the numerical evolution of the field according to Eq. (j3J) are shown in Fig. [2j 
The left panel shows the field evolution along <fi and ip directions when A = 10 _10 mpi, / = mpi, 
and the initial misalignment 9 = 7r/10000. Note that the initial misalignment is very small 
from t/> = 0, the field rolls towards a saddle point and hence <fi settles down at a minimum first. 
This is clearly seen from several trajectories shown in the right panel of Fig. [2J Fig. [3] illustrates 
the dependence of the number of e-folds on the energy scale A and the decay constant /. 
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Figure 2: (Left) evolution along <j) and ijj directions when A = 10 _10 mpi, / = m-pi, and 9 = 
7r/10000, and (right) several field trajectories for which the parameters are given by the same 
as the left panel (red), A = 10 -9 mpi, / = 2 x 10 _1 mpi, and 9 = 7r/100 (green, magnified by 3 
times), and A = 10 -8 mpi, / = 5 x 10 -2 mpi, and 9 = tt/10 (blue, magnified by 6 times). 

One interesting thing is that locked inflationary phase [7] hardly occurs. This could be seen 
from the fact that there is no interaction term in the potential, e.g., such as \<p 2 ip 2 . Thus the 
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Figure 3: Number of e-folds depending on A and /, with 9 = 7r/100. 



oscillation of, say, is not transmitted to the effective mass of ifj to hold the field at a saddle 
point. Note that even if we include an interaction term of the form 



this term is suppressed by m P1 and does not lead to any significant change. Also we note that 
the dependence of the number of e-folds on the initial misalignment 9 [8] is very weak. We are 
considering / < mj>\, and the slow-roll parameter \rj\ > 0(1), see Eq. ([3]). Hence the field is very 
quickly rolling, and this makes the total number of e-folds more or less the same irrespective 
of the initial misalignment 9. In Table El we present several results depending on 9. 



Table 1: Several results depending on different initial misalignment 9. Parameters are given by 
A = 10~ 10 mpi and / = m-p\. Note that the number of e-folds varies less than 1. 

Now let us turn our attention to the perturbations produced during inflation. At the early 
stage of inflation when the field is very close to a maximum, we can split the field into the 



4 In general, the field follows curved and possibly complicated trajectory on the potential. If the number of 
e-folds is large enough by e.g. some mechanism discussed in Section [4] and the perturbations produced when 
the field is around a maximum is far outside the horizon, then this decomposition may complicate the matter. 
In this case, one can just write [9] 




(9) 



9 tt/100 tt/40 tt/20 tt/10 tt/8 tt/5 
N 30.4465 30.2241 30.2100 30.1144 30.3783 30.0687 




(p = a/0 2 + -02 , 



(10) 




to calculate the power spectrum of the adiabatic perturbations after inflation. 
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and 6 as Eq. Note that near a maximum, U(l) symmetry is good enough. The radial 
component (p thus plays the role of the "inflaton" , yielding adiabatic perturbations, while the 
angular component 9 is orthogonal to the field trajectory and its fluctuations correspond to 
the isocurvature perturbations. Therefore the radial fluctuations follow the form of those from 
the usual inverted quadratic potential pU], but the angular ones are very close to those of 
massless scalar field and the spectrum is nearly scale invariant [El [H] . Now, we can follow the 
5N formalism [9] to write the power spectrum of the primordial density perturbations as 

2 

- 1 v. . 



dN\ 2 fdN\ 2 , . 

where it is clear that both the adiabatic and the isocurvature components contribute to the final 
perturbations. For the isocurvature perturbations to dominate the primordial perturbations, 
we require that 

w) >> W) ■ (12) 

Since the adiabatic component is boosted by a factor of generally we need (dN/d8) 2 ^> 1, 

i.e., the number of e-folds is highly dependent on where the field is rolling towards. But as 
we have seen in Table [31 the number of e-folds is more or less the same irrespective of the 
misalignment. Thus (dN/d$) 2 < 1 so the contribution of the isocurvature component is 
negligible, and the density perturbation spectrum is completely dominated by the contribution 
of the inflaton component which is highly scale dependent. 

Now let us consider the case when ^ fy, i.e. the potential is not symmetric but 
"squeezed". Then, the field tends to roll along the squeezed direction first irrespective of 
the initial misalignment, because that direction is (much) steeper than the other one. Hence 
the isocurvature perturbations, i.e. the perturbations orthogonal to the field trajectory, become 
irrelevant and the "inflaton" perturbations are dominant. Let us assume that ip direction is 
squeezed so that = afy with a > 1. With the rotated slow- roll parameters 

riuu = VH c °s 2 ® + sin2 & 

Vss = sin2 & + V-ipip c °s 2 & > ( 13 ) 

where a and s denote the inflaton and the isocurvature components respectively, the evolution 
of the isocurvature perturbations is given bj|f| [12] 



S ~ (-2e + rjcra - ri ss )HS 

HS 



-2 e + ^-^^|i(cos 2 0-sin^) 

« 2 fl 



a 2 — 1 m|j 



-2-—^]^, (14) 



5 Note that although this equation is valid under slow-roll approximation, nevertheless it still gives physically 
interesting insights. 



5 



which is always negative. Intuitively, as a direction becomes more and more squeezed, the 
potential becomes closer to the single field case where only the inflaton perturbations exist. 



4 Remedies 

In the previous section, we have seen that the number of e-folds obtained during inflation 
driven by two decoupled string axion fields is generally not enough to solve many cosmological 
problems. Also the power spectrum of density perturbations is completely dominated by the 
radial contribution which is not scale invariant. The naive inflation model of the two string 
axions is thus not sufficient and we need some amendments to improve the situation. In this 
section we briefly summarise some known remedies. 

One obvious modification is to add more contribution to H so that the field receives greater 
friction and can roll down the potential slowly. One would be tempted to add a positive 
cosmological constant, but this may cause two difficulties. When the cosmological constant is 
too large, it may render the universe to expand forever even after the field has settled down at 
a minimum. Also, since the observed dark energy density is of O (lO" 12 eV 4 ) , the magnitude of 
the cosmological constant one would like to add is severely constrained. It is possible, however, 
to increase H by adding more fields and make them rolling slowly for many e-folds [13]. Then 
Eq. (j4j) becomes 

« +3 *<-%t« + w^mr =0 - (15) 

where the contribution H by other fields is dependent on models: for example, if we adopt 
N-flatiorj§ [TJ] and add k additional axions, we have 



3m P1 # 2 = V 



ml<Pl, (16) 

where we assume that each axion field is displaced not too far from the minimum. It is possible 
that by some of the added fields further perturbations are generated after inflation [16], but 
the parameter space is heavily restricted for this mechanism to work properly |17j . 

Another clear way we can take is to make the axion decay constant / larger than mpi and 
follow the original wisdom of natural inflation [18] . One difficulty we meet is that because it is 
required that / ^> mpi to satisfy the slow-roll conditions, / is outside the region of the effective 
field theory description and is out of our control. An interesting alternative to build a model 
where the above problem could be evaded is to identify the extra component A$ of a gauge field 
in a 5D theory with the inflaton [19] . Because this model is based on higher dimensional theory, 
the effective decay constant / e g = {T.Tig^R)^ 1 ^> m P i, where R is the radius of the circle where 
the extra dimension is compactified, is reliable. Another way of constructing / > mp\ is to 
incorporate two different gauge groups, making the effective decay constant larger than mpi by 
imposing a symmetric relation between the couplings of the axions to the gauge groups [20J. 



6 In a subsequent study |15j . we found that although the detail depends on the underlying mass distribution, 
roughly 1/10 of the total fields contribute at the end of the inflationary phase. 
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In Fig. H] several results of e-folds depending on / and 9 are shown. Note that when the field 
rolls towards a saddle point, generally it acquires extra e-folds but the scale invariance of the 
perturbation spectrum is broken when the field leaves the saddle point [21 J. To avoid this, we 
should put strong constraints on the parameters [6]. 
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Figure 4: Plot of N versus / and 9. The energy scale A is set to be 10 _4 mpi. We can see that 
as / becomes bigger, we obtain larger iV as expected. Also note that the dependence on 9 is 
not significant. 

Also note that when / ^> mpi, the minima are separated by a distance far larger than 
mpi. Hence, the extrema which intervene between those minima become the sources of eternal 
topological inflation [22] . For example, a saddle point which separates two minima along (say) 
ijj direction, it becomes a domain wall and in its vicinity we have an eternally inflating regime. 

One fundamental question we can ask at this point is whether it is possible to find regimes 
of moduli space where / ^> mp\ in string theory. It seems that, unfortunately, there exists no 
consistent scenario yet [23J although not all the possibilities are explored. Building a consistent 
scenario where a large decay constant could arise is still an interesting and challenging topic. 

5 Summary 

We have investigated a simple inflation model which consists of two non-interacting axions. 
The number of e-folds A" achieved while the field rolls off from a maximum is dependent on 
both the energy scale A and the decay constant /: we obtain larger A" with smaller A and 
with bigger /, but generally A" is not sufficient to solve cosmo logical problems. Interestingly, in 
spite of the rich geometry of the potential, the initial misalignment 9 hardly affects A" because 
of the absence of the interaction term, which also leads to no locked inflationary phase. The 
perturbation spectrum is highly scale dependent, because the scale invariant spectrum which 
arises from the isocurvature component is greatly suppressed. To alleviate this situation, we 
may either increase H by adding more axion fields, or make / bigger than mp\ which is still 
unclear in the context of string theory. 
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